Holographic dark energy models: A comparison from the latest observational data 

Miao Li, 1 2 Q Xiao-Dong Li, 3 ' Shuang Wang, 4 - and Xin Zhang 5 - 2 J3 

'institute of Theoretical Physics, Chinese Academy of Science, Beijing 100080, China 
2 Kavli Institute for Theoretical Physics China, 
Chinese Academy of Sciences, Beijing 100080, China 
3 Interdisciplinary Center for Theoretical Study, 
University of Science and Technology of China, Hefei 230026, China 
4 Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China 

5 Department of Physics, College of Sciences, 
Northeastern University, Shenyang 110004, China 

The holographic principle of quantum gravity theory has been applied to the dark energy (DE) 
problem, and so far three holographic DE models have been proposed: the original holographic 
dark energy (HDE) model, the agegraphic dark energy (ADE) model, and the holographic Ricci 
dark energy (RDE) model. In this work, we perform the best-fit analysis on these three models, by 
using the latest observational data including the Union+CFA3 sample of 397 Type la supernovae 
(SNIa), the shift parameter of the cosmic microwave background (CMB) given by the five-year 
Wilkinson Microwave Anisotropy Probe (WMAP5) observations, and the baryon acoustic oscillation 
(BAO) measurement from the Sloan Digital Sky Survey (SDSS). The analysis shows that for HDE, 
y 2 . = 465.912; for RDE, v 2 . = 483.130; for ADE, y 2 . = 481.694. Among these models, HDE 

I mm ' ' A mm ' ' A mm & ' 

model can give the smallest x m ,„■ Besides, we also use the Bayesian evidence (BE) as a model 
selection criterion to make a comparison. It is found that for HDE, ADE, and RDE, A In BE = -0.86, 
-5.17, and -8. 14, respectively. So, it seems that the HDE model is more favored by the observational 
data. 

PACS numbers: 98.80.-k, 95.36,+x 



I. INTRODUCTION 



Observations of Type la supernovae (SNIa) [1], cosmic microwave background (CMB) [2] and large 



scale structure (LSS) 



all indicate the existence of mysterious dark energy (DE) driving the current accel- 



erating expansion of the universe. The most obvious theoretical candidate of dark energy is the cosmological 



'Electronic address: mli@itp.ac.cn 
^Electronic address: renzhe@mail.ustc.edu.cn 
^Electronic address: swang@mail.ustc.edu.cn 
^Electronic address: zhangxin@mail.neu.edu.cn 



2 



constant A, which can fit the observations well, but is plagued with the fine-tuning problem and the coinci- 
dence problem [4]. Numerous other dynamical DE models have also been proposed in the literature, such 
as quintessence Ja], phantom |fj], ^-essence |7], tachyon [8], hessence [sj], Chaplygin gas flo| . Yang-Mills 
condensate [11], ect. 

I I 

Actually, the DE problem may be in essence an issue of quantum gravity 111211 . However, by far, a 
complete theory of quantum gravity has not been established, so it seems that we have to consider the effects 
of gravity in some effective field theory in which some fundamental principles of quantum gravity should 
be taken into account. It is commonly believed that the holographic principle H 1 3M is just a fundamental 



principle of quantum gravity. Based on the effective quantum field theory, Cohen et al. il4fl pointed out that 
the quantum zero-point energy of a system with size L should not exceed the mass of a black hole with the 
same size, i.e. L}p vac < LM 2 ^, where p vac is the quantum zero-point energy density, and Mpi = 1 / VSnG 
is the reduced Planck mass. This observation relates the ultraviolet (UV) cutoff of a system to its infrared 
(IR) cutoff. When we take the whole universe into account, the vacuum energy related to this holographic 
principle can be viewed as dark energy (its energy density is denoted as hereafter). The largest IR cutoff 
L is chosen by saturating the inequality, so that we get the holographic dark energy density 



Pde 



3c 2 M 2 pl L- 2 , 



(1) 



where c is a numerical constant characterizing all of the uncertainties of the theory, and its value can only be 
determined by observations. If we take L as the size of the current universe, for instance the Hubble radius 



H 1 , then the dark energy density will be close to the observational result. However, Hsu H 15f1 pointed out 



this yields a wrong equation of state (EOS) for DE. 

Li lllql suggested to choose the future event horizon of the universe as the IR cutoff. This is the original 



holographic dark energy model. Subsequently, Cai proposed 111711 that the age of the universe can be chosen 



as the IR cutoff, and this model is called agegraphic dark energy model. A new version of this model 



replacing the age of the universe by the conformal age of the universe was also discussed 11811 . in order to 



avoid some internal inconsistencies in the original model. Furthermore, Gao et al. [19] proposed to consider 
the average radius of the Ricci scalar curvature as the IR cutoff, and this model is called holographic Ricci 
dark energy model. (For convenience, hereafter we will call them HDE, ADE, and RDE, respectively.) 
Although these three models have been studied widely 
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22], so far no comparison of them has been 



made. It would be very interesting to constrain the holographic dark energy models by using the latest 
observational data, and then make a comparison for them by using the proper model selection criterion. 
This will be done in this work. 

This paper is organized as follows: In Section 2, we briefly review the holographic DE models. In 
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Section 3, we present the method of data analysis, as well as the model selection criterion. In Section 4, 
we show the results of the cosmological constraints, and we also use the Bayesian evidence (BE) to make 
a comparison. Section 5 is a short summary. In this work, we assume today's scale factor ao = 1, so the 
redshift z satisfies z = a~ l - 1; the subscript "0" always indicates the present value of the corresponding 
quantity, and the unit with c = h = 1 is used. 



II. MODELS 



In this section, we shall briefly review the holographic DE models. For a spatially flat (the assumption 
of flatness is motivated by the inflation scenario) Friedmann-Robertson- Walker (FRW) universe with matter 
component p m and dark energy component pj e , the Friedmann equation reads 



or equivalently, 



3M 2 PI H 2 =p m + Pde , (2) 



Hq \ 1 - Side 



where H = a /a is the Hubble parameter, £l m o is the present fractional matter density, and £1^ = = 3! ^t H 2 
is the fractional dark energy density. In this work, for simplicity, we shall not discuss the issue of interaction 
between dark matter and dark energy, so we have 

p n + 3H Pm = , (4) 



Pde + 3H(l +w de )p de = 0, (5) 
where the over dot denotes the derivative with respect to the cosmic time t, and Wd e is the EOS of DE. 



A. The HDE model 



For this model, the IR cutoff is chosen as the future event horizon of the universe, 

da' 



r°° dt' r° 

L — a I — - a I 

Jt a J a 



Ha 



(6) 



Taking derivative for Eq. (Q} with respect to x = In a and making use of Eq. (O, we get 

, dpde V^rfT 



(7) 



4 



Combining Eqs. © and ([7]), we obtain the EOS for HDE, 



1 2 , 

Wde = —Z ~ T" S&de- (8) 
3 3c 



Directly taking derivative for Qde = c 2 /(H 2 L 2 ), and using Eq. (©, we get 

J&de 



n' d =m de [e-l + ^-^\, (9) 



where e = -HIE 1 = -H'/H. From Eqs. ©, ©, ©, and ©, we have 

3„ „ , 3 n de n 3 J 2 



e = ^(1 + w^O*) = 1 - ^ - (10) 
2 2 2 c 



for this case. So, we have the equation of motion, a differential equation, for Q de , 



_q;,, = n,,,n -q (/< )(i + ^V^|- (ii) 



Since ^ = -(1 + z)^, we get 



" Q "' V -(\ +z y 1 a de (l-Q de )h + -^ de \. (12) 



Solving numerically Eq . (fT2l) and substituting the corresponding results into Eq . (0, the key function E(z) 
can be obtained. It should be mentioned that there are two model parameters, il m o and c, in the HDE model. 

B. The ADE model 

Since there are some internal inconsistencies in the original ADE model, we will discuss the new version 
of ADE model which suggests to choose the conformal age of the universe 

/dt C da „„ s 

7 = Jm- (13) 

as the IR cutoff, so the density of ADE is 

p de = 3n 2 M 2 Pl rj- 2 . (14) 

To distinguish from the HDE model, a new constant parameter n is used to replace the old parameter c. 
Taking derivative for Eq. (fT4l with respect to x and making use of Eq. (TT3T ). we get 

p' de = ~2p de ^L. (15) 
na 



This means that the EOS of ADE is 



2 'sj^lde 

Wde = ~l + - ■ (16) 

in a 
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Taking derivative for D.de - n 2 /(H 2 r] 2 ), and considering Eq. (TT3T ). we obtain 



Q' =2Q d A e -^^\. (17) 
\ na } 

In this case, we have 

3 3 3 3/2 

e = -(1 + vv^O^) = - - -H de + (18) 
2 2 2 «a 

Hence, we get the equation of motion for Q. c / e , 

2 V^rfe 



a' de = o«fc(i - n&) 3 — — , (19) 

1 « a 



and this equation can be rewritten as 

dO-de 



-Q de (l - Q de ) 13(1 + z)" 1 - 2 - Vn^J 



(20) 



n. 

As in 11811 . we choose the initial condition, Qdeiztni) = n 2 (l + zmT 2 74, at zm - 2000, then Eq. (1201 can be 
numerically solved. Substituting the results of Eq. (1201) into Eq . ©, the key function E(z) can be obtained. 
Notice that once n is given, Q m Q = 1 - Qde(z = 0) can be natural obtained by solving Eq.(l20l. so the ADE 
model is a single-parameter model. 

C. The RDE model 

For a spatially flat FRW universe, the Ricci scalar is 

K = -6(H + 2H 2 ). (21) 

n 

As suggested by Gao et al. 11911 . the energy density of DE is 



3a 

Pde 



(H + 2H 2 ) = —ft, (22) 



8^G V > \6nG 

where a is a positive numerical constant to be determined by observations. Comparing to Eq. £[]), it is seen 
that if we identify the IR cutoff L with -ft/6, we have a = c 2 . As pointed out by Cai et al. [23], the RDE can 
be viewed as originated from taking the causal connection scale as the IR cutoff in the holographic setting. 
Now the Friedmann equation can be written as 

H 2 - — p m0 e- 3x + a I - — + 2H 2 , (23) 

and this equation can be further rewritten as 

E 2 = a m0 e-^ + al\^ + 2E 2 \, (24) 
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where E = H/Hq. Solving this equation, and using the initial condition Eq = E(to) = 1, we have 

E(z) = ( + z) 3 + (1 - + zf- 2 A V2 . (25) 

\2 — a 2 - a ) 

There are also two model parameters, Q m Q and a, in RDE model. 



III. METHODOLOGY 
A. Data analysis 

In the following, we constrain the model parameters of these three DE models by using the latest ob- 
servational data including the Union+CFA3 sample of 397 SNIa, the shift parameter of the CMB given 
by the five-year Wilkinson Microwave Anisotropy Probe (WMAP5) observations, and the baryon acoustic 
oscillation (BAO) measurement from the Sloan Digital Sky Survey (SDSS). 

First, we consider the latest 397 Union+CFA3 SNIa data, the distance modulus Hobsizd, compiled in 
Table 1 of 12411 . This dataset add the CFA3 sample to the 307 Union sample 12511 . The theoretical distance 
modulus is defined as 



Hthizt) - 5 log 10 D L {zi) + juo, (26) 
where hq = 42.38 - 5 log 10 h with h the Hubble constant Hq in units of 100 km/s/Mpc, and 

is the Hubble-free luminosity distance H^di (here di is the physical luminosity distance) in a spatially flat 
FRW universe, and here 6 denotes the model parameters. The^ 2 for the SNIa data is 

397 

2 ,o\ - V WobsiZi) ~ HthiZi)] 



^ _ y \MobsW) ~ fJMZM 



(=1 i 



where Hobsizd and cr,- are the observed value and the corresponding lcr error of distance modulus for each 
supernova, respectively. The parameter fio is a nuisance parameter but it is independent of the data and the 
dataset. Following [26], the minimization with respect to //o can be made trivially by expanding the^- of 
Eq. (l28l with respect to /uo as 

X 2 SN (9) = A{9) - 2 m B(6) + n 2 C, (29) 



where 



A(6) = f ^MzM^d^l3 , ( 30) 
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B{9) = £ 



Hobsizd- Hth(zi\ Ho = 0,8) 



(31) 



i=l 



tr- 



ial 



Evidently, Eq. (1281 has a minimum for fiQ = B/C at 



(32) 



2 B(6) 2 
f SN (0) = A{0) - --f 



(33) 



Since X^ Nmin = X^N m i n > instead minimizing x^ N one can minimize x 2 SN which is independent of the nui- 
sance parameter fiQ. 

Next, we consider the constraints from CMB and LSS observations. For the CMB data, we use the CMB 
shift parameter R; for the LSS data, we use the parameter A of the BAO measurement. It is widely believed 
that both R and A are nearly model-independent and contain essential information of the full CMB and LSS 

|3Q 



BAO data. The shift parameter R is given by 

Jo 



dz' 

E(z')' 



(34) 



where the redshift of recombination z re c - 1090 [29]. The shift parameter R relates the angular diameter 
distance to the last scattering surface, the comoving size of the sound horizon at z rt > r and the angular scale of 

281. The measured value of 



the first acoustic peak in CMB power spectrum of temperature fluctuations 127 
R has been updated to be R obs = 1.710 + 0.019 from WMAP5 The parameter A from the measurement 
of the BAO peak in the distribution of SDSS luminous red galaxies is defined as 13011 

-i2/3 



a s titter" 3 



Zb Jo 



dz' 



E(z') 



where z b = 0.35. The SDSS BAO measurement 1301] gives A obs = 0.469 («. s /0.98) 



n-0.35 



(35) 



0.017, where the 



scalar spectral index is taken to be n s = 0.960 as measured by WMAP5 12911 . Notice that both 7? and A are 
independent of Ho, thus these quantities can provide robust constraint as complement to SNIa data on the 
holographic dark energy models. The total x 1 is given by 



2 _ ~2 2 2 

X - XsN + XcMB + XBAO 



(36) 



where x\ N is given by Eq. d33l , and the latter two terms are defined as 



2 (R ~ Robs) 

XcMB - 



(37) 



cr 



R 



8 



and 



2 _ (A A b s ) 2 
XBAO - 2 ' ( --'°- , 



where the corresponding lcr errors are cr R = 0.019 and cta = 0.017, respectively. As usual, assuming the 
measurement errors are Gaussian, the likelihood function 

L oc e~* 2/2 . (39) 

The model parameters yielding a minimal^ 2 and a maximal X will be favored by the observations. 



B. Model Comparison 

For comparing different models, a statistical variable must be chosen. The;^ 2 . is the simplest one which 
is widely used. However, for models with different number of parameters, the comparison using x 1 may not 
be fair, as one would expect that a model with more parameters tends to have a lower x 1 - In this work, we 
will use the BE as a model selection criterion. The BE of a model M takes the form 

BE = J £(d\6, M)p(6\M)dG, (40) 

where -C(d\9, M) is the likelihood function given by the model M and parameters 9 (d denotes the data), and 
p(6\M) is the priors of parameters. It is clear that a model favored by the observations should give a large BE. 
Since BE is the average of the likelihood of a model over its prior of the parameter space and automatically 
includes the penalties of the number of parameters and data, it is more reasonable and unambiguous than 
the^. in model selection 13 111 . The logarithm of BE can be used as a guide for model comparison [31], 
and we choose the ACDM as the reference model: A In BE = lnBE,„ 0£ / e / - InBE^coM- 



IV. RESULTS 

We will show the results in this section. In FigfJJ we plot the contours of 68.3% and 95.4% confidence 
levels in the Q m o - c plane, for the HDE model. The best-fit model parameters are Q m o = 0.277 and 
c = 0.818, corresponding to x 2 m in = 465.912. For 68.3% confidence level, Q m0 = 0.277+J}^ and c _ 
0.818+°-^; for 95.4% confidence level, n m0 - 0.277^7, and c = 0.818+° J^. 

In Fig 12 we plot the x 1 an d the corresponding likelihood X of the ADE model. The best-fit model 
parameter is n = 2.807, giving x 2 min = 481.694. For 68.3% confidence level, n = 2.807+^^; for 95.4% 
confidence level, n = 2.807+°-™. 
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FIG. 1: Probability contours at 68.3% and 95.4% confidence levels in Q. m o - c plane, for the HDE model. The round 
point denotes the best-fit values, O m o = 0.277 and c = 0.818, corresponding toxt- = 465.912. 
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FIG. 2: The^ 2 and the corresponding likelihood X. of the ADE model. 



In Figj3j we plot the contours of 68.3% and 95.4% confidence levels in the Q m o - a plane, for the RDE 



model. The best-fit model parameters are Q m o - 0.324 and a = 0.371, corresponding to^ 



2 _ 



483.130. 



For 68.3% confidence level, Q m0 = 0.324+°-°^, and a = 0.371+°°|; for 95.4% confidence level, Q m0 = 
324+ 040 and a = 371 +0.037 

-0,036' -0.038- 

Therefore, among these three holographic DE models, the HDE model can give the smallest x 2 . 



As 



a comparison, we also fit the ACDM model to the same observational data, and find that the minimal 
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FIG. 3: Probability contours at 68.3% and 95.4% confidence levels in Q„,o - a plane, for the RDE model. The round 
point denotes the best-fit values, Q,„o = 0.324 and a = 0.371, corresponding to k;. = 483.130. 




FIG. 4: Comparison of the observed H(z), as square dots, with the predictions from the holographic DE models. 



X min = 467.775 for the best-fit parameter O m o = 0.274. Note that both HDE model and RDE model are 
two-parameter model, while ADE model is a single-parameter model. So, it is not suitable to choose the 



X 2 min as a model selection criterion. In Table|H we list the results of A In BE for three holographic DE models. 



It is seen that although HDE model performs a little poorer than ACDM model, it performs better than ADE 
model and RDE model. Therefore, among these three holographic DE models, HDE is more favored by the 
observational data. 

In FigJU we further compare the observed expansion rate H(z) 13211 with that predicted by these three 
holographic DE models (for each model, the best-fit values determined by SNIa+CMB+BAO analysis are 
taken). Notice that the area surrounded by two dashed lines shows the 68% confidence interval [2], and a 



11 



TABLE I: The results of A In BE for three holographic DE models. 



Model 


HDE 


ADE 


RDE 


A In BE 


-0.86 


-5.17 


-8.14 



DE model would be disfavored by the observation if it gives a curve of H(z) that falls outside this area. It 
is seen that among these three holographic DE models, only the curve of H(z) predicted by the HDE model 
falls inside this confidence interval. This result verifies our conclusion, from another perspective. 

V. SUMMARY 

In this work, we perform the best-fit analysis on three holographic DE models, by using the latest 
observational data including the Union+CFA3 sample of 397 SNIa, the shift parameter of the CMB given 
by the WMAP5 observations, and the BAO measurement from the SDSS. The analysis shows that for HDE, 
y 2 . = 465.912; for RDE, v 2 . = 483.130; for ADE, y 2 . = 481.694. Among these three DE models, 

A. mm ' ' n. mm ' ' a. ffufi o ' 

the HDE model can give the smallest v 2 . , which is even smaller than that given by the ACDM model. 
Moreover, we use the BE as a model selection criterion to make a comparison. Although HDE model 
performs a little poorer than ACDM model, it performs better than ADE model and RDE model. Therefore, 
among these three DE models, HDE is more favored by the observational data. Finally, adopting the best-fit 
values determined by SNIa+CMB+BAO analysis, we compare the observed expansion rate H(z) with that 
predicted by these three holographic DE models. It is found that only the curve of H(z) predicted by HDE 
model can fall inside the 68% confidence interval. This verifies our conclusion from another perspective. 
It would be interesting to give tighter constrains on DE models by adding other cosmological observa- 



tions, such as gamma-ray bursts 



330, Chandra x-ray observation [34], and some old high-redshift objects 



1 35J]. This will be studied in a future work. 
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